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Abstract 



JlLi' We use the stochastic quantization method to study systems with complex valued path integral 



weights. We assume a Langevin equation with a memory kernel and Einstein's relations with 
colored noise. The equilibrium solution of this non-Markovian Langevin equation is analyzed. 
We show that for a large class of elliptic non-Hermitian operators acting on scalar functions on 
Euclidean space, which define different models in quantum field theory, converges to an equilibrium 

00 \ state in the asymptotic limit of the Markov parameter r — > oo. Moreover, as we expected, we 

cn . obtain the Schwinger functions of the theory. 
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1 Introduction 



The program of stochastic quantization, first proposed by Parisi and Wu [1] , and the stochastic 
regularization was carried out for generic fields defined in fiat, Euchdean manifolds. A brief 
introduction to stochastic quantization can be found in Refs. [2] [3] [4], and a complete review 
is given in Ref. [5]. Recently Menezes and Svaiter [6] implemented the stochastic quantization 
in the theory of self-interacting scalar fields in a static Riemannian manifold and also a manifold 
with a event horizon, namely, the Einstein and the Rindler manifold. First, these authors solved a 
Langevin equation for the mode coefficients of the field, then they exhibit the two-point function 
at the one-loop level. It was shown that it diverges and to regularize the theory they used a 
covariant stochastic regularization. The presence of the Markov parameter as an extra dimension 
allows the authors to implement a regularization scheme, which preserves all the symmetries of 
the theory under study. It is clear that the stochastic quantization program can be implemented 
without problems, if it is possible to perform the Wick rotation, obtaining a real Euclidean action. 

The picture that emerges from the above discussion is that the implementation of the stochastic 
quantization in curved background is related to the following fact. For static manifold, it is 
possible to perform a Wick rotation, i.e., analytically extend the pseudo- Riemannian manifold 
to the Riemannian domain without problem. Nevertheless, for non-static curved manifolds we 
have to extend the formalism beyond the Euclidean signature, i.e., to formulate the stochastic 
quantization in pseudo-Riemannian manifold, not in the Riemannian space (as in the original 
Euclidean space) as was originally formulated. See for example the discussion presented by Huffel 
and Rumpf [7] and Gozzi [8]. In the first of these papers the authors proposed a modification 
of the original Parisi- Wu scheme, introducing a complex drift term in the Langevin equation, to 
implement the stochastic quantization in Minkowski spacetime. Gozzi studied the spectrum of the 
non-self-adjoint Fokker-Planck Hamiltonian to justify this program. See also the papers [9] [10]. 
Of course, these situations are special cases of ordinary Euclidean formulation for systems with 
complex actions. 

The main difference between the implementation of the stochastic quantization in Minkowski 
spacetime and in Euclidean space is the fact that in the latter case the approach to the equilibrium 
state is a stationary solution of the Focker-Planck equation. In the Minkowski formulation, the 
Hamiltonian is non-Hermitian and the eigenvalues of such Hamiltonian are in general complex. The 
real part of such eigenvalues are important to the asymptotic behavior at large Markov time, and 
the approach to the equilibrium is achieved only if we can show its positive semi-definiteness. The 
crucial question is: what happens if the Langevin equation describes diffusion around complex 
action? Some authors claim that it is possible to obtain meaningful results out of Langevin 
equation diffusion processes around complex action. Parisi [11] and Klauder and Peterson [12] 
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investigated the complex Langevin equation, where some numerical simulations in one-dimensional 
systems was presented. See also the papers [13] [14]. We would also like to mention the approach 
developed by Okamoto et al. [15] where the role of the kernel in the complex Langevin equation 
was studied. 

We would like to remark that there are many examples where Euclidean action is complex. 
The simplest case is the stochastic quantization in Minkowski spacetime, as we discussed. Other 
situations are systems, as for example QCD, with non- vanishing chemical potential at finite tem- 
perature; for SU{N) theories with N > 2, the fermion determinant becomes complex and also the 
effective action. Complex terms can also appear in the Langevin equation for fermions, but a suit- 
able kernel can circumvent this problem [16] [17] [18]. Another case that deserves our attention is 
the stochastic quantization of topological field theories. One of the peculiar features within these 
kind of theories is the appearance of a factor i in front of the topological actions in Euclidean space. 
In these topological theories, the path integral measure weighing remain to be e*'^, even after the 
Wick rotation. An attempt to use a Markovian Langevin equation with a white noise to quantize 
the theory fails since the Langevin equation will not tend to any equilibrium at large Markov 
parameter. In the literature there are different proposes to solve the above mentioned problem. In 
a pure topological Chern-Simons theory, Ferrari et al. [19] introduced a non-trivial kernel in the 
Langevin equation. Other approach was developed by Menezes and Svaiter [20]. These authors 
proved that, using a non-Markovian Langevin equation with a colored random noise, it is possible 
to obtain convergence towards equilibrium even with an imaginary Chern-Simons coefficient. An 
interesting application of this method can be found on Ref. [21], where a Langevin equation with 
a memory kernel was introduced in order to obtain the Schwinger functions for the self-interacting 
scalar model. In conclusion, although several alternative methods have been proposed to deal 
with interesting physical systems where the Euclidean action is complex [22] [23] [24] [25], these 
methods do not suggest any general way of solving the particular difficulties that arise in each sit- 
uation. Here, we wish to report progress in the stochastic quantization of theories with imaginary 
action, introducing a memory kernel. 

It is the purpose of the present paper to use the method of the stochastic quantization to 
study systems with complex valued path integral weights. We assume a Langevin equation with a 
memory kernel and Einstein's relations with colored noise [26]. We show that for a large class of 
elliptic non-Hermitean operators which define different models in quantum field theory converges 
in the asymptotic limit of the Markov parameter r — *> oo, and we obtain the Schwinger functions 
of the theory. In section II, we briefiy discuss the Parisi-Wu stochastic quantization for the case 
of free scalar field. In section III we implement the stochastic quantization for scalar field with 
complex action using a non-Markovian Langevin equation. Conclusions are given in section IV. 
In this paper we use h = c = ks = I. 



2 Stochastic quantization for the free scalar field theory: 
the Euclidean case 

In this section, we give a brief survey of the stochastic quantization. This technique in flat space- 
time with trivial topology can be summarized by the following steps. First, starting from a fleld 
defined in Minkowski spacetime, after analytic continuation to imaginary time, the Euclidean 
counterpart, i.e., the field defined in an Euclidean space, is obtained. Second, it is introduced 
a monotonically crescent Markov parameter, called in the literature "fictitious time" and also a 
random noise field ri{T,x), which simulates the coupling between the classical system and a heat 
reservoir. It is assumed that the fields defined at the beginning in a d-dimensional Euclidean space 
also depends on the Markov parameter, therefore the field and a random noise field are defined in a 
((i+ l)-dimensional manifold. One starts with the system out of equilibrium at an arbitrary initial 
state. It is then forced into equilibrium assuming that its evolution is governed by a Markovian 
Langevin equation with a white random noise field [27] [28] [29]. In fact, this evolution is de- 
scribed by a process which is stationary, Gaussian and Markovian. Finally, the n-point correlation 
functions of the theory in the {d + l)-dimensional space are defined by performing averages over 
the random noise field with a Gaussian distribution, that is, performing the stochastic averages 
{(p{Ti,Xi)Lp{T2,X2)--.'^{Tn,Xn) ),,• The n-poiut Schwiuger functions of the Euclidean d-dimensional 
theory are obtained evaluating these n-point stochastic averages ( (/?(ri,Xi)</?(r2,X2)...¥?(T„, a;„) )^ 
when the Markov parameter goes to infinity (r -^ oo), and the equilibrium is reached. This can 
be proved in different ways for the particular case of Euclidean scalar field theory. One can use, 
for instance, the Fokker-Planck equation [30] [31] associated with the equations describing the 
stochastic dynamic of the system. A diagrammatical technique [32] has also been used to prove 
such equivalence. 

After this brief digression, let us consider a free neutral scalar field. The Euclidean action that 
usually describes such theory is 

So[^] = Jd'^x(^^{d^r + lml^'{x)y (1) 

The simplest starting point of the stochastic quantization to obtain the Euclidean field theory 
is a Markovian Langevin equation. Assume a fiat Euclidean (i-dimensional manifold, where we are 
choosing periodic boundary conditions for a scalar field and also a random noise. In other words, 
they are defined in a d-torus Q = T^. To implement the stochastic quantization we supplement the 
scalar field if{x) and the random noise 'r]{x) with an extra coordinate r, the Markov parameter, 
such that if{x) — > ip{T,x) and 'r]{x) -^ ?7(r, x). Therefore, the fields and the random noise are 
defined in a domain: T'^ x R^^\ Let us consider that this dynamical system is out of equilibrium. 



being described by the following equation of evolution: 



d_ 



V9(r,x) 



5^0 



6ip{x] 



\lf(x) = Lp{T,x) + 'HV^I •^/l 



(2) 



where r is a Markov parameter, 77 (r, x) is a random noise field and Sq is the usual free Euclidean 
action defined in Eq. (1). For a free scalar field, the Langevin equation reads 



^V'(r,x) 



-A + ml)ip{T,x) +r]{T,x), 



(3) 



where A is the rf-dimensional Laplace operator. The Eq. (3) describes a Ornstein-Uhlenbeck 
process and we are assuming the Einstein relations, that is: 



(?7(r,x))^ = 0, 
and for the two-point correlation function associated with the random noise field 

( r]{r, x) r]{r', x') ), = 2S{r - r') S\x - x'), 



(4) 



(5) 



where { ■■■)n means stochastic averages. The above equation defines a delta-correlated random 
process. In a generic way, the stochastic average for any functional of (p given by F[ip] is defined 

by 



{F[^])v 



J[dri]F[ip]exp 



■\jd'^xJdTr]^{T,x) 



I [dv] exp 



— J J d'^x j dT7f{T,x) 



(6) 



Let us define the retarded Green function for the diffusion problem that we call G{t — r',x — x'). 
The retarded Green function satisfies G{t — r', x — x') = if r — r' < and also 



d_ 



+ {-A^ + ml 



G{t -t',x- x') = 6\x - x')6{r - r'). 



(7) 



Using the retarded Green function and the initial condition ip{T,x)\r=o = 0, the solution for Eq. 
(3) reads 

^(r, x) = r dr' f d'^x' G{t -t\x- x')r]{T\ x') . (8) 

Let us define the Fourier transforms for the field and the noise given by (/^(r, A;) and ?7(r, A;). We 
have respectively 



V^(t, k) 



{2n\ 



dxe (p{T,x), 



(9) 



and 

viT,k) = -^ f d''xe-^'^v{T,x). (10) 

[ZTTj 2 J 

Substituting Eq. (9) in Eq. (1), the free action for the scalar field in the ((i + l)-dimensional space 
writing in terms of the Fourier coefficients reads 

So[v{k)]\ip{k)=^{T,k) = 1^ (fk(f{T,k){k^ + ml)ip{r,k). (11) 

Substituting Eq. (9) and Eq. (10) in Eq. (3) we have that each Fourier coefficient satisfies a 
Langevin equation given by 

— V9(r, k) = -(P + ml)^{r, k) + r^{r, k). (12) 

In the Langevin equation the particle is subject to a fluctuating force (representing a stochastic 
environment), where its average properties are presumed to be known and also the friction force. 
Note that the "friction coefficient" in the Eq. (12) is given by (A;^ + rn^). 
The solution for Eq. (12) reads 

V9(r, k) = exp {-{^ + ml)T) V9(0, A;) + j' dr'exp (-(A;' + rnD^r - r')) r]{r', k). (13) 

Using the Einstein relation, we get that the Fourier coefficients for the random noise satisfies 

{Vir,k)), = (14) 

and 

( rjiT, k)r]iT', k') ), = 2 5{t- t') 5\k + k'). (15) 

It is possible to show that ( (/?(r, k)(f{T', k') ),y|r=r' = D{k, k'; r, r') is given by: 

Dik; T, r) = i2n)''S'ik + k')—^-^(l - exp (-2r(A;2 + m^))) . (16) 

where we assume r = r'. Therefore, for r ^ oo we recover the Euclidean two-point function. 

The self-interacting theory is beyond the scope of this paper, however it can be carried out in a 
straightforward way. In the next section we present a modification of the Langevin equation that 
allows us to treat systems with complex Euclidean actions. Although non-trivial, it is intuitively 
obvious that the method can be extended to interacting field theory. 



3 Stochastic quantization for complex actions 

As an application of the ideas discussed previously, in this section we show how it is possible to 
quantize a theory with a complex action using a non-Markovian Langevin equation. We will be 
following similar steps as in Ref. [21]. Consider the following Euclidean action: 

S = ]^j d'^x^K^, (17) 

with the following Markovian Langevin equation: 

d 

—(p{T,x) = -K'f + r]{T,x). (18) 

where K is an elliptic operator (with some minor changes, our proof in this section can be made 
to hyperbolic operators). The function y? is a scalar field, for simplicity, but we can generalize our 
results to fields of higher spin. If we let K to be non- Hermit ian, the action in Eq. (17) becomes 
complex. There are many approaches in the literature to deal with complex actions; one of them 
is to employ a complex Langevin equation, separating the field in a real part and in a imaginary 
part, Re{ip) = ipi and Im{ip) = ip2 [5]. With this approach, we get two Langevin equations, one 
for each of the two fields ipi and (f2- Another one is to use a modified Langevin equation: 



Q^<f{T,x) = - J d yK{x,y) ^ . J y(y)=y(r,;;) + ??(r, x), (19) 

where a subsequent change in the second moment of the noise field is 

( ri{T, x) riir', x'))r, = 26{t - r) k{x, x') . (20) 

With these modifications, we may choose an appropriate kernel: 

k{x,x') = Kl6{x-x'), (21) 

so the Langevin equation becomes 

d 
— V9(r, x) = -K.K^ifi + r]{T, x). (22) 

We see that we get a "bosonized" version of the Langevin equation given by Eq. (18) and the 
problem with the convergence towards an equilibrium disappears, since K.K'^ is a Hermitian 
operator. This is the prescription usually employed in the literature to deal with the stochastic 
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quantization of fermions [16] [17] [18]. In fact, the root of this problem hes in the fact that there 
exists no classical analogue of fermion fields. 

Another approach can be used as well. Let us consider the following non-Markovian Langevin 
equation: 

d r 5 S 

-—ip{T,x) = - dsM^{T - s) A ^(x)=v{s,x)+v{t,x), (23) 

OT JO (p{X) 

where M\ is a memory kernel and the stochastic random field ?](r, x) satisfies the modified Ein- 
stein's relations 

( vir, x) r]ir', x') ), = 2MA(|r - r'\) S\x - x'). (24) 

In this case where Mx{\t — r'\) has a width in the fictitious time, the description is Gaussian in 
spite of being non-Markovian. For the case of Euclidean free scalar field theory we have that the 
generalized Langevin equation reads 

d r 

—-(p(T,x) = - dsM/^ir - s)Kip(s,x) +r](T,x). (25) 

OT JO 

We shall prove in this section that this method leads to convergence towards equilibrium, even 
though we have a complex Langevin equation. 

We can introduce a mode decomposition such as 



ip{T,x) = / dfi{n)ipn{T)un{x) (26) 

and 

ri{T,x) = / dfi{n)r]n{T)un{x), (27) 

where the measure fi{k) depends on the metric we are interested in. Each Fourier coefficient ipn 
obeys a (non-Markovian) Langevin equation given by 

d r 

--V9„(r) = -A„ / rfsMA(r-s)v9„(s) + ?7„(r), (28) 

OT JO 

where A„ is an eigenvalue of the operator K and r7„(r) obeys 

( Vn{r) r]^iT') )^ = 2MA(|r - r'|) 5\n, n'). (29) 

Following Fox [33] [34], we define the Laplace transform of the memory kernel: 

Miz) = / dTMA(T)e-'\ (30) 

JO 



With the initial condition ipn{T~)\T=o = 0, the solution of the Eq. (28) reads: 

(^„(r)= / dT'G^{r-r')7]^{T'), (31) 

Jo 

where using the step function 0{t), the Green function Gn{T) is defined by: 

Gr.{T)=n,,{r)e{T). (32) 

The fi„(r) function that appears in Eq. (32) is defined through its Laplace transform: 

^n{r) = [z + A„M(z)) . (33) 

From Eq. (31) and the modified Einstein relations, we get that the free scalar correlation function 
Dn{T,r') is given by: 

Dnir,T') = 

= 2S'^{n,n') ds ds' Gn{T - s)Gn{T' - s') M^{\ s - s' \) 
Jo Jo 

= 2S'^{n,n') r ds f rfs'fi„(r - s) fi„(r' - s') ^^a(| s - s' |). (34) 



JO 



To proceed we have to write Dn{T, r') in a simplified way. Note that the double Laplace transform 
of the right hand side is given by: 

/•OO /"OO rr rr 

/ dre-'^ dr'e"^ ds ds'f]„(r - s) f]„(r' - s') Ma(| s - s' |) = 
Jo Jo Jo Jo 

= nn{z)^n{z') ds ds'e-''' e-''M^{\s-s' \). (35) 

JO Jo 

Now, with simple manipulations, we get: 

d.rd.'e-'--'e-"M^(|.-. -I) - '"'' + ':"''' ■ (36) 

Jo z + z' 

Therefore, we get the identity: 

/•OO roo rr rr' 

/ dre-'^ dr'e-'^ ds rfs'f]„(r - s) f]„(r' - s') Ma(| s - s' |) = 
Jo Jo Jo Jo 

^n„Wn„(.')( ^^'-~'t"f'-'' )^ (37) 

V z + z' I 



Remembering Eq. (33), we can show that: 

\ Z + Z J ^n \ Z + Z J 

So, in parallel with result given by Eq. (36), we finally obtain a very simple expression for Dn{T, r') 
in terms of r2„(r). We have 

D„(r, r') = y 5''in, n') (^fi„(| r - r' |) - f]„(r) f]„(r')) . (39) 

Now, we need an expression for the memory kernel in order to investigate the convergence of Eq. 
(39). A series of kernels were proposed in the literature: 

Mjr\r) = ^A^(A^ I r |)™exp(-A2 | r |). (40) 

For simplicity, we shall take the case for m = 0. Then, from Eq. (30), Eq. (33) and Eq. (40), and 
applying the inverse Laplace transform [35] , we obtain the following expression for the fi-function: 

^n{r) = (f «i^h(^) +cosh(^)) ^^p(-^4")' ^^^^ 

where we have defined a quantity /9 given by: 

f3 = x + ty, (42) 

where 

X = ,/ZIiiIH, (43) 



2(A4 + a„ + |;|)' 
and, finally: 



y = ^nhiT^. m^^ (44) 



'{K^ + anY + al (45) 

with a = — 2A^A„ and we have written a as a = aji + iai, where aji and ai are real quantities. 
Similarly, we have, for the Green function: 

Gn{r) = (f si^h(^) +cosh(^)') exp(-r^)^(r). (46) 



Now, in order to have convergence, we must demand that Gn -^ when r — » oo. Using hyperbohc 
identities we see that, in order to have G„ — ;► when r — *> oo, we shall have: 






(Al 



/\2 



2A2 ' 



(47) 



where we have written the eigenvalues as A„ = A^ + iA^. Since A is, in principle, arbitrary, 
we see that the condition given by Eq. (47) does not seem to pose any serious restrictions on 
the eigenvalues of the operator K. However, we have another restriction. Note that, with this 
prescription, the function M{x — y; r) defined in the Ref. [21] as 

M{x -y,T)= r ds MAir - s) G{t - s, X - y), (48) 

Jo 

where G{t — s, x) is the retarded Green function for the diffusion problem, whose Fourier transform 
is given by: 



A2 

M{k,T) = — 



2 9A4 + /32 

. f 3A2 
4 exp(^ — r 



8A^ 



2A^cosh^ + 



Pr 



3A4 



— I sinh ( — 
2/ V 2 



will no longer be Hermitian. But, if we impose the following restriction: 

lA^I < 5A2, 



(49) 



(50) 



its real part will remain positive, which implies that the real part of the eingenvalues of the Fokker- 
Planck Hamiltonian defined therein are positive, assuring, therefore, convergence to equilibrium, 
i.e., the system reaches its ground state. It seems that Eq. (47) and Eq. (50) can be imposed 
simultaneously without restrictions on A. In particular, we are allowed to take arbitrarily large 
values for A, which would imply that our approach works for almost any value for the eigenvalues 
A„. It can be easily proved that the presence of the zero modes for massless theories destroys the 
convergence to equilibrium. 

From the results above, it is easy to see that the free two-point function is given by: 



An 



— sinh ( — 

V 2 



cosh ( — 
V 2 



r — T 



+ cosh 



Pir - t') 



A2 
expi-- 



[T — T 



—— sinh ( I + cosh ( 

/9 V 2 y V 2 



exp 



A^ 



r + T 



(51) 
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For r = r', we get: 

Dnir,T) = ^5\n,n') U - [^ sinh(^) + cosh(^)') expl-A^)") . (52) 

So, in the limit r — ;► oo, we obtain the following result: 

D^{T,T) = ^5\n,n'). (53) 

Therefore, apart from an unimportant constant, we have obtained that, in the asymptotic 
limit r — s> oo, we have reached convergence towards the expected equilibrium, and the two-point 
Schwinger function was obtained. 



4 Conclusions and perspectives 

There are many examples where the Euclidean field theory is defined for an imaginary action. 
Since in this case the path integral weight is not positive definite, the stochastic quantization in 
this situation is problematic. Parisi and Klauder proposed complex Langevin equations [11] [12], 
and some problems of this approach are the following. First of all, complex Langevin simulations 
do not converge to a stationary distribution in many situations. Besides, if it does, it may converge 
to many different stationary distributions. The complex Langevin equation also appears when the 
original method proposed by Parisi and Wu is extended to include theories with fermions [16] [17] 
[18]. The first question that appears in this context is if make sense the Brownian problem with 
anticommutating numbers. It can be shown that, for massless fermionic fields, there will not be a 
convergence factor after integrating the Markovian Langevin equation. Therefore the equilibrium 
is not reached. One way of avoiding this problem is to introduce a kernel in the Langevin equation 
describing the evolution of two Grassmannian fields. 

In this paper, we have used the method of the stochastic quantization to study systems with 
complex valued path integral weights. We assumed a Langevin equation with a memory kernel 
and Einstein's relation with colored noise. The equilibrium solution of such Langevin equation 
was analyzed. We have shown that for a large class of elliptic non-Hermitean operators which 
define different models in quantum field theory converges in the asymptotic limit of the Markov 
parameter r -^ oo, and we have obtained the Schwinger functions of the theory. Although non- 
trivial, the method proposed can be extended to interacting field theory with complex actions. 
The generalization of the method for this situation is under investigation by the authors. 
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